Probing the Electronic Structure of Bilayer Graphene by Raman Scattering 
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The electronic structure of bilayer graphene is investigated from a resonant Raman study using 
different laser excitation energies. The values of the parameters of the Slonczewski-Weiss-McClure 
model for graphite are measured experimentally and some of them differ significantly from those 
reported previously for graphite, specially that associated with the difference of the effective mass 
of electrons and holes. The splitting of the two TO phonon branches in bilayer graphene is also 
obtained from the experimental data. Our results have implications for bilayer graphene electronic 
devices. 
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After the development of the phenomenological 
Slonczewski-Weiss-McClure (SWM) model for the elec- 
tronic structure of 3D graphite 50 years ago Jj, i2,] , many 
experimental works were devoted to evaluate the param- 
eters of this standard model [1, 3- In particular, the 
weak interaction between two layers in graphite, ascribed 
to van der Waals forces, can be described by the in- 
terlayer SWM parameters, that give, for example, the 
difference between the effective masses of electrons and 
holes near the Dirac point (K point in the graphite Bril- 
louin zone) This work shows that, by performing 
Raman scattering experiments in bilayer graphene with 
many different laser excitation energies, we can probe 
its electronic structure by obtaining experimental values 
for the SWM parameters. Our values for some inter- 
layer parameters differ significantly from values measured 
previously for graphite 0, 3], showing that open issues 
about graphite can now be experimentally revisited using 
graphene samples Furthermore, while the unbiased 
bilayer graphene is a zero-gap semiconductor, a biased 
bilayer is a tunable gap semiconductor by electric field 
effect [1, 0] ■ Hence, the development of bilayer graphene- 
based bulk devices depends on the detailed understand- 
ing of its electronic properties. 

Fig. [1] shows the atomic structure of a bilayer graphene. 




FIG. f : (a) Atomic structure of bilayer graphene. The A 
atoms of the two layers are over each other, whereas the B 
atoms of the two layers are displaced with respect to each 
other. The SWM constants 70, 71, 73 and 74 label the corre- 
sponding pair of atoms associated with the hopping processes, 
(b) First Brillouin zone of monolayer graphene, showing the 
high symmetry points F, K, K' and M. 



in which we can distinguish the two nonequivalent atoms 
A and B in each plane giving rise to a unit cell with four 
atoms. Since this unit cell is the same for graphite in the 
Bernal stacking structure, we can describe the electronic 
spectrum of bilayer graphene in terms of the SWM model 
for graphite [H, by determining the parameters 70, 71, 
73 and 74, that are associated with overlap and trans- 
fer integrals calculated for nearest neighbors atoms. The 
pair of atoms associated with these parameters are indi- 
cated in the atomic structure of a bilayer graphene shown 
in Fig. [IJa). These parameters, that are fundamental for 
the electronic processes in the system, are only roughly 
known to this date. 

The graphene samples used is this experiment were ob- 
tained by a micro-mechanical cleavage of graphite on the 
surface of a Si sample with a 300 nm layer of Si02 on the 
top tSl]- The bilayer flakes were identified by the slightly 
color change from monolayer graphene in an optical mi- 
croscope, followed by a Raman spectroscopy character- 
ization using the procedure described by Ferrari et al. 
[3]. For the Raman measurements, we used a Dilor XY 
triple monochromator in the back scattering configura- 
tion. The spot size of the laser was ~1 /xm using a 100 
X objective and the laser power was kept at 1.2 mW in 
order to avoid sample heating. Raman spectra were ob- 
tained for 11 different laser lines of Ar-Kr and dye lasers 
in the range 1.91-2.71 eV. 

Recently, Ferrari et al. P showed that Raman spec- 
troscopy can be used to identify the number of layers in a 
graphene sample and, in particular, to clearly distinguish 
a monolayer from a bilayer graphene sample. Figure [2] 
shows the Raman spectra of the monolayer [Fig. [Ifa)] 
and bilayer [Fig. [D^b)] graphene samples, where the most 
prominent features are the G and G' Raman bands 0]. 
The G' band of the monolayer graphene is nicely fitted by 
just one Lorentzian, whereas four Lorentzian peaks are 
necessary to fit the G' band of the bilayer graphene, in 
agreement with the previous Raman studies of graphene 
systems [l[iS[ill. 

The Raman spectra of both the monolayer and bi- 
layer graphene have been measured with many different 
laser energies in the visible range. Figure [3] shows the 
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FIG. 2: (a) Raman spectrum of the monolayer graphene and 
(b) Raman spectrum of the bilayer graphene performed with 
the 2.41 eV laser line. The G' band of bilayer graphene is fit 
with four different Lorentzians with FWHM of ~ 24 cm~^. 
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FIG. 3: (a) Laser energy dependence of the G'-band energy 
for a monolayer graphene. (b) Laser energy dependence of 
the positions of the four peaks that comprise the G'-band of 
bilayer graphene. These four Raman peaks originate from the 
Pii, Pi2, P21 and P22 DR processes illustrated in Fig. |4l The 
black squares are the experimental data and the full lines are 
the fitted curves (see discussion in the text). 



laser energy dependence of the G'-band frequency for the 
monolayer sample [Fig. [Slja)] and for each one of the four 
peaks that comprise the G' band for bilayer graphene 
[Fig.Hb)]. 

The origin of the G' band in all graphitic materials 
is due to an inter-valley double-resonance (DR) Raman 
process involving electronic states near two nonequiva- 
lent K points (K and K') in the 1** Brillouin zone of 
graphene [see Fig . [Hb)], and two phonons of the iTO 
branch [i, d m,. As a result of the angular depen- 
dence of the electron-phonon matrix elements and 
the existence of interference effects in the Raman cross- 
section [l^ , the main contribution for the G'-band comes 
from the particular DR process that occurs along the F- 
K-M-K'-F direction, in which the wavevectors k and fc' 
of the two intermediate electronic states in the conduc- 
tion band (measured from the K and K' points, respec- 
tively) are along the KF and K'F directions, respectively 
Q. Therefore, the wavevector q of the phonons involved 
in this specific process is along the KM direction and 
is related to the electron wavevectors by the condition 

In the case of the bilayer graphene, where the two 
graphene layers are stacked in a Bernal configuration Q , 
the TT-electrons dispersion in the valence and in the con- 
duction band splits into two parabolic branches near the 
K point [iBl, as shown in Fig. H) In this figure, the upper 
and lower branches of the valence band are labeled as tti 
and 7r2, respectively. The lower and upper branches of 
the conduction band are called ■k\ and ■k^^ respectively. 
Along the high symmetry F-K-M-K'-F direction, these 
branches belong to different irreducible representations 
of the P63/mmc space group and, therefore, only the 
7ri<=^7r* and 7r2«=^7r2 optical transitions between the va- 
lence and conduction bands are allowed. Therefore, there 
are four possible intervalley DR processes involving elec- 
trons along the F-K-M-K'-F direction that lead to the 
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FIG. 4: Schematic view of the electron dispersion of bilayer 
graphene near the K and K' points showing both tti and H2 
bands. The four DR processes are indicated: (a) process Pn, 
(b) process P22, (c) process P12, and (d) process P21. The 
wavevectors of the electrons (fci, k2, k'l and k'2) involved in 
each of these four DR processes are also indicated. 



observation of the four peaks of the G'-band of bilayer 
graphene 

The four DR processes are represented in Fig. 01 In 
process Pn [Fig. (DJa)], an electron with wavevector fci is 
resonantly excited from the valence band tti to the con- 
duction band tt^ by absorbing a photon with energy ■ 
This electron is then resonantly scattered to a state with 
wavevector k[ by emitting a phonon with momentum qn 
and energy E^^ . Finally, the electron is scattered back to 
state ki by emitting a second phonon, and it recombines 
with a hole producing a scattered photon with an energy 
Es — El — 2£'p^. The phonon wavevector qn, measured 
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from the K point and along the KM direction, is given 
by qii = ki + k[. 

Fig. \Mi>) shows the DR process P22, which involves 
an electron that is optically excited between the n2 and 
TTj branches. The energy of the associated phonon is 
Ep"^ and its wavevector is given by q22 = fc2 + ^2- Fig- 
ures Sl^c) andSJd) show processes P12 and P21, that in- 
volve electrons with wavevectors fci, ^2 and A:2, k'l, re- 
spectively, which belong to different electron branches. 
The wavevectors of the phonons associated with pro- 
cesses P12 and P21 are given by qi2 = ki + k'2 and q2i 
— k2 + k[, respectively. 

The two iTO phonon branches of a bilayer graphene 
along the PKM line belong to the irreducible represen- 
tations Ti and T3. The scattering of an electron in 
the conduction band between states around K and K' 
has to satisfy the electron-phonon selection rule, so that 
the allowed transitions are ttI^ttI or tt2^'^2 foi' the Ti 
phonon and ■K2'^t^i for the T3 phonon. Since processes 
Pii and P22 involves electrons states around K and K' 
which belong to the same electronic branch, the associ- 
ated phonons belong to the Ti phonon branch. On the 
other hand, phonons involved in processes P12 and P21 
belong to the T3 branch. 

From Fig. |4] one can see that the phonon associated 
with the Pii process has the largest wavevector (qn). 
Since the energy of the iTO phonon increases with in- 
creasing g, the highest frequency component of the G'- 
band of a bilayer graphene is related to the Pn process. 
On the other hand, P22 process is associated with the 
smallest phonon wavevector q22, and gives rise to the 
lowest frequency component of the G'-band of a bilayer 
graphene. The two intermediate peaks of the G'-band 
are associated with processes P12 and P21. 

In order to analyze the experimental results shown 
in Fig. mb), we must find a relation between the elec- 
tronic and the phonon dispersions of a bilayer graphene, 
according to the DR Raman process. The electronic 
dispersion of the bilayer graphene can be described in 
terms of the standard SWM model for graphite (see 
e.g. Eq. (2.1) in Ref. [l|) and using the full tight- 
binding dispersion introduced originally by Wallace [l7| . 
Along the K — Y direction this amounts to replacing 
cr by 7o[2cos(27r/3 - fca\/3/2) + l] in McClure's ex- 
pressions. Here k is measured from the K point and 
a = 1.42 A is the in plane nearest neighbor carbon dis- 
tance. Since there is no fc^-dependence we may set F = 1 
and 72 = 75 = (72 and 75 correspond to 3'"'' layer inter- 
actions in graphite). We have verified that the parameter 
A does not make any noticeable difference in our results 
and will be ignored in our analysis. With these simplifi- 
cations, the bands in the bilayer are parameterized by the 
four parameters 70, 71, 73 and 74. A more detailed ac- 
count on this approach to the band structure calculation 
of the graphene bilayer can be found in Ref. [2]| . 

In fact, along the high symmetry K — T direction, the 
4x4 matrix factorizes and the dispersion of the four 



bands are given by: 

- (-71 - «3^T - e+)/2, (la) 

E^, = (7i+i'3T-e-)/2, (lb) 

E.i = (-7i-i;3^ + e+)/2, (Ic) 

E^'^ = (7i+«3a + e-)/2, (Id) 

where 

e± = v/(7i-^'3^)'+4(l±«4)V2, (2) 

and Vj = 7j/7o. 

In order to obtain the dependence of the phonon energy 
E'^J on the laser energy El, let us consider a generic 
process P^j, where i, j = 1 or 2, which describes the 
four processes shown in Fig. |4l In the first step of this 
process (electron- hole creation), the incident photon is 
in resonance with the electronic states in the valence and 
conduction bands at the ki point. Thus, the laser energy 
El can be written as: 

EL = E^*{h)-E^^{ki). (3) 

Eq. [3] allows us to determine the momentum ki of the 
electron excited in the process. 

The electron is then resonantly scattered from the 
vicinity of the K point to the vicinity of the K' point 
by emitting a (iTO) phonon with energy E^^ given by: 

E^ {h + k'^)^ E^^ ih)^ m ) (4) 

Assuming that we know the iTO phonon dispersion near 
the K point [A + B{ki + kj)] as well as the bands involved 
(Eqs. 1 and [2), Eq. [?] uniquely determine the momen- 
tum of the scattered electron k'j. We then compute the 
phonon energy E^^ , that is directly related to the Raman 
shift for this specific Pij process, obtained with a given 
laser energy El- Finally, we perform a least-squares fit to 
determine the parameters 70 , 71 , 73 and 74 of the model 
(Eqs. 1 and [5]) that give the best fit of the dispersion of 
the four G' peaks in bilayer graphene shown in Fig. [31 
The four solid lines in Fig. ^h) represent the best fitting 
of the experimental E^^ versus El data. 

We have also tried to fit, unsuccessfully, our experi- 
mental data for the four DR processes taking only 70, 71 
and 73. Therefore, in order to get a good fitting, the pa- 
rameter 74 value, that is associated with the splitting of 
the two G' intermediate peaks, has to be included. For 
the dispersion of the iTO phonon branches near the K 
point, we could not fit satisfactorily the dispersion data 
in Fig. [3]Jb) considering the same phonon dispersion for 
the four Py processes. The best fit was obtained when we 
considered different dispersions for the two iTO phonon 
branches of bilayer graphene. Table [J shows the param- 
eters obtained for each phonon branch, which exhibit a 
linear dispersion near the K point. Notice that the differ- 
ence between these two phonon branches corresponds to 
a splitting of about 3 cm~^, which is in close agreement 
to that reported previously [1, |Tl[ . 
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TABLE I: Values obtained for the phonon dispersion of the 
two iTO phonon branches of bilayer graphene near the K 
point, where E^—Au+'Biiq corresponds to the Pn and P22 
processes and Ep^ =Aij+'Bijq to the P12 and P21 processes. 









Bij 


153.7 (meV) 
1238 (cm'i) 


38.5 (meVA) 
310 (cm-l A) 


154.0 (meV) 
1241 (cm-i) 


38.8 (meVA) 
313 (cm-l A) 



Table |TT] shows the 7 values obtained experimentally. 
The parameter 70, associated with the in-plane nearest- 
neighbor hopping energy, is ten times larger than 71, 
which is associated with atoms from different layers along 
the vertical direction (see Fig. [T|). These values are in 
good agreement with the previous angle resolved photoe- 
mission spectroscopy (ARPES) measurements in bilayer 
graphene [13] ■ The resolution of our experiment allows, 
however, the measurement of weaker hopping parameters 
(73 and 74) , that are beyond the current resolution of 
ARPES. The value of 71 is about three times larger than 
73 and 74, both associated with the interlayer hopping 
not along the vertical direction (see Fig. [T]) . 

The corresponding parameters found experimentally 
for graphite are also shown in Table HIl [3, ij] . The param- 
eters 7o and 71 for bilayer graphene are slightly smaller 
than those for graphite. This difference is more accentu- 
ated for the parameter 73. Notice that our value for 73 
is in good agreement with that reported theoretically by 
Tatar and Rabii IT§\. On the other hand, our value of 74 
for bilayer graphene is significantly higher than the value 
for graphite measured by Mendez et al. in a magnetore- 
flection experiment. [l9|. Moreover, this parameter is 
specially important since it is related to the difference 
of electron and hole effective masses in the valence and 
conduction bands This is particularly important in 



the context of the biased bilayer graphene We must 
stress that our value for 74 is in good agreement with the 
calculations performed by Tatar and Rabii [11] and by 
CharHer et al. 



TABLE IL Experimental SWM parameters (in eV) for the 
band structure of bilayer graphene. The parameters for 
graphite are taken from Refs. 0- 





70 


71 


72 


73 


74 


75 


76 = A 


Bilayer graphene 


2.9 


0.30 


n/a 


0.10 


0.12 


n/a 


n/a 


Graphite 


3.16 


0.39 


-0.02 


0.315 


0.044 


0.038 


0.008 



In summary, from the resonant Raman study of the 
G'-band of bilayer graphene using several laser excita- 
tion energies, we have been able to probe the dispersion 
of electrons and phonons of this material near the Dirac 
point. From the fitting of the experimental data using the 
SWM model, we have obtained experimental values for 
hopping parameters for bilayer graphene. The parame- 
ters 73 and 74 found here differ significantly from the val- 
ues previously measured for graphite [1, [j, [l^ , showing 
that open issues about graphite can now be revisited us- 
ing graphene samples. The slight difference between the 
experimental and theoretical data in Fig. [3)^b), specially 
for lower laser energies, mi ght be ascribed to many-body 
effects in bilayer graphene [23j . Finally, in a future work 
we will study DR processes involving holes as well as the 
biased bilayer graphene. 
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